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The diffuse Galactic 7-ray data from the region of the Galactic Center has been collected by the 
LAT instrument on the Fermi Gamma- Ray Space Telescope. In this paper we argue that it may be 
able to provide an unambiguous evidence of originating, in addition to known astrophysical sources, 
^>«^ ' from dark matter annihilations in the halo, independently of the mass and other properties of the 

I dark matter particle. We also show that the recently released high precision data from mid-latitudes 

, is already providing an upper bound, albeit still a weak one, on the cuspiness of the dark matter 

' density profile as a function of the mass of the dark matter assumed to be a stable neutralino of 

^ minimal supersymmetry. 
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I. INTRODUCTION 



The evidence for the existence of dark matter (DM) in the Universe is strong and mounting but its nature remains 
' unknown. It is generally believed that it is made up of some unknown weakly interacting massive particle (WIMP) 
for which there are a whole host of possible hypothetical candidates predicted in various extensions of the Standard 
Model of particle physics, with the lightest neutralino of supersymmetry (SUSY) being perhaps the most popular 
one 0]. Nevertheless, it is clear that the mass and other properties of DM can only be established experimentally. A 
1^, whole range of detection strategies have been developed to search for DM. A number of underground detectors are 
P ' currently in operation trying to directly detect a faint signal from a Galactic halo WIMP scattering off a target and, if 
successful, they would provide perhaps the most unambiguous signal of DM. Indirect searches look for exotic products 
^ • of WIMP annihilation (or decay, if unstable) in the Sun's or Earth's core, in the Galactic halo or the Galactic Center 

■ (GC) [H. Finally, WIMP particles are Ukely to be produced at the LHC. 
Following last year's launch of the Fermi Gamma-Ray Space Telescope (Fermi), 7-ray window on the sky has 

' received a major boost with new unprecedented quality data on pulsars gamma-ray-bursters 0], diffuse radiation, 
'f' , etc. The diffuse radiation in particular is of much interest since, under favorable conditions, it may reveal a measurable 

■ contribution from annihilation products of dark matter [l|. Such favorable conditions may exist in the region of the 
ly-j , GC where DM density is believed to be enhanced, and Fermi data from that region are eagerly awaited by the dark 
^-H ' matter community. In the meantime, high precision data from Fermi LAT for diffuse emission with photon energies 

of 0.1 GeV to 10 GeV from the region of intermediate Galactic latitudes 10° < \b\ < 20° and the full range of the 
Galactic longitudes {0 < I < 360°) has recently been presented 0, 
' In this paper we make two points. The first is that, by comparing upcoming Fermi data on 7-ray flux at different 
, angular distances from the Galactic Center, one may be able to unambiguously infer its DM origin, independently of 
^ ' the dark mater mass and other properties. The second point is that, the recently released Fermi data for mid-latitude 
. ^ , 7~rays already provides a constraint on the profile of the Galactic DM halo as a function of DM mass assumed to be 

■ the lightest neutralino of supersymmetry. 

[ The paper is organized as follows. In Sec. we briefly review the formalism used to compute 7-rays from DM 
. - - . annihilation and list some popular DM halo profiles used subsequently in Sec. IIIII to suggest tests of DM signatures 
in the GC and to derive, in Sec. HVl an upper bound on the cuspiness of the DM halo profile. In Sec. |V] we briefly 
summarize our results. 



II. DIFFUSE GAMMA-RAYS FROM GALACTIC DARK MATTER ANNIHILATION 

The differential diffuse 7-ray flux originating from WIMP pair-annihilation, and arriving from a direction at an 
angle ip from the GC is given by [l| 

where m^^ denotes the WIMP mass and a^v is a product of the WIMP pair-annihilation cross section into a flnal state 
i times the pair's relative velocity and dN1^/dE^ is the differential 7-ray spectrum (including a branching ratio into 
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photons) following from the state i. Here we consider the total contribution from both the continuum, resulting from 
cascade decays of all kinematically allowed final state SM fermions and combinations of gauge and Higgs bosons, and 
from photon lines coming from one loop direct neutralino annihilation into 77 and jZ. The integral over the square 
of the DM mass density is taken along the line of sight (l.o.s.) I' at an angle ip between us and the GC. It is related 
to the Galactic longitude (1) and the Galactic latitude (6) by cos ip = cos I cos b. 

It is convenient to separate factors depending on particle physics and on halo properties by introducing the dimcn- 
sionless quantity 



.5kpc VO.3 GeV/cm- 



l.o.s. 



The flux arriving from the angle ^ is further averaged over the solid angle AO representing point spread function 
angle of the detector or some wider angle of the sky, depending on the situation, and one defines the quantity 

iJW)An^^J^J^'J(^^')- (3) 

The total flux from a solid angle AQ centered on the angle tp and integrated over photon energy from an energy 
threshold i?th is then given by 



$.^(A0) = / dE-yd^-y/dEj{E-y,An). (4) 

One of the crucial ingredients in the analysis is the distribution of dark matter in and around the Milky Way (MW) . 
Unfortunately, despite much effort and progress, this remains rather poorly known, especially in the inner part of the 
MW. For this reason, in this analysis we will consider some popular halo profiles with a strongly varying cuspiness at 
small Galactic radius r. 

Several popular profiles can be parametrized by Q 

[l + (ro/r,)"]^ 

Px (r) = Po — , (5) 

(r/ro)^[l + (r/r,n^ 

where po — 0.3 GeV/cm'^ represents the local DM density. For example, in the Navarro, Frenk and White (NFW) 
model To = 8.0 kpc, = 20.0 kpc, a = 1, /3 = 3 and 7 = 1 Another model that has recently become favored by 
numerical simulations of large structures is the Einasto model [1, described by 

2 f°' — f°' 

p(r)=poexp[-^-^], (6) 

where in this case po = 0.066 GeV/cm'^ and = 25 kpc, while the index a takes the range j < a < ^ [T^, [ll[, 
with the best-fit value of a = 0.17 [l2]. In this paper, we will consider the best-fit case (called thereafter "Einasto"), 
although one should remember that increasing a within the above range gives a whole range of profiles with decreasing 
cuspiness. The Burkert model is an example of a realistic profile with a basically flat DM density distribution in 
the inner region of the MW. It is parametrized by 

, , [1 + jro/an)] [1 + (ro/a„)^] 
[1 + (r/a„)J [1 + {r/any\ 

where po = 0.34 GeV/cm^ and a„ = 11.68 kpc. As the last profile, we consider the model of Klypin, et al., [l3| as 
providing the most divergent profile in the GC. It is based on the NFW model but it is fitted to the data from the 
MW and takes into account the effect of angular momentum exchange between baryons and dark matter, and in this 
sense may be considered as more applicable to our Galaxy than other profiles. The inner radius density profile for 
the Klypin, et al., model is ~ r^^-^ which is steeper than the Einasto and NFW profiles. Close to the solar radius all 
the models become quite similar. As the Klypin, et al., model results from the fits to the data, no analytic expression 
is available. 

The above DM density profiles are presented in Fig. [TJ We will use them as illustrative examples as they represent 
a large variation in cuspiness at small r, which will be perhaps the single most important property in determining 
DM contribution to diffuse 7-radiation. 



3 




FIG. 1: Some popular dark matter halo profiles as a function of radius r. 



III. SIGNATURES OF DARK MATTER IN 7-RAY SPECTRUM FROM THE GALACTIC CENTER 



Ideally, one of the most convincing signatures of the annihilating DM origin of the diffuse 7-radiation that Fermi's 
data could produce, would be the flux with a spectrum that falls off with the angle tp from the GC. This is because 
it is proportional to p^, with the DM density in most models decreasing with r, or -0, compare Eq. ([1]) and Fig. [1] 
Since the flux itself does depend on the DM particle mass and its other unknown properties, as well as the photon 
energy E-y, we propose instead to consider the ratio 

which follows from Eqs. (Il])-([3|) and which is clearly only dependent on the dark matter mass density distribution in 
the halo. The square dependence of along the l.o.s. to a large degree cancels out, but the one in the transverse 
direction does not, and provides a genuine effect of DM annihilation in the region of the GC. The ratio ([5]) is shown 
in Fig. [2] for three representative halo models and for two different angular resolutions typical for Fermi LAT. (In the 
case of the more cuspy profiles, in order to avoid a divergent behavior, we set a cutoff radius of ~ 10^'^kpc.) 

Observing the ratio R'^^ shown in Fig. [5] would constitute a clear and rather unambiguous signal of the DM 
origin of the diffuse 7-radiation. Indeed, it is highly unlikely that any other, known or unknown, source of diffuse 
7-radiation would have a distribution producing a similar fall-off. At least we are not aware of any.^ The astrophysical 
sources of the diffuse flux in the photon energy range of interest are relatively well understood. The main contribution 
is caused by primary cosmic ray (CR) interactions with interstellar hydrogen and helium atoms, producing tt^'s which 
in turn decay via 7r° 77 with a peak at ~ 70 MeV, with a symmetric distribution. Other sources include: 

inverse Compton fluxes and bremsstrahlung due to CR scatterings of electrons and positrons off the interstellar 
medium, isotropic background (e.g., extragalactic radiation) and point sources. These contributions are well modelled 
with GALPROP which in general reproduces current data remarkably well. 

Of course, as we said above, finding evidence for DM in the measurement of the ratio ^ would require a highly 
cuspy profile at small r and also relatively small astrophysical background. While the former is poorly known and 



^ We stress that the behavior presented in in Fig. [2] applies to stable DM, whose pair annihilation is proportional to p|. The flux from 
decaying DM instead depends on the DM density linearly and would produce a much flatter spectrum which would be probably harder 
to distantangle from astrophysical contributions. 
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FIG. 2: The ratio -R^s^/dis^ = {^(11^0))'' (^^e GC) versus V for AH = 10~^sr (dashed) lO'^sr (sohd) for the halo models 
considered in the paper. The ratio does not depend on the DM particle mass nor its other microscopic properties. For the 
Galactic latitude b — the angle coincides with the Galactic longitude, = ±Z, while for I — ^ coincides with tp — ±6 (up 
to 90°). 



remains a possibility, the latter is known to be hugely important and is likely to overshadow, or even make it too 
difficult, to observe the behavior shown in Fig. [21 Firstly, the DM is expected to be strongly enhanced only within the 
rather small range of some 2 — 3° around "0 = 0. Secondly, there are many known (and also likely many unresolved) 
point sources in the direction of the GC. However, the DM signal should be basically spherical while the astrophysical 
component is expected to be disk-like.^ For this reason it will be of crucial importance to examine the ratio shown 
in Fig. [5] in the special cases of the Galactic plane {b — 0) and the vertical plane / = 0, where the astrophysical 
component should show a different angular behavior but the DM part should remain the same. The test is certainly 
going to be very challenging but should be attempted as it would provide perhaps a single most convincing test of 
the DM origin of the measured ratio. Even if it is possible to measure the ratio Eq. ([5]) only for some ranges of ip, 
this already could provide some vital information pointing towards the DM origin of the effect. 

Furthermore, if the measured ratio shows the behavior shown in Fig. [2l one could attempt to actually infer DM 
profile, at least to some degree. In particular, if the ratio drops rapidly at small ip^ this would imply a very cuspy 
profile, which of course will also produce a much larger absolute flux at small r, and for this reason will be much 
easier to detect. 

Therefore, a measurement of the ratio Eq. ([8]), if it confirms the behavior shown in Fig. [21 will provide a convincing 
signature of existence of DM in the region of the GC, independently of the DM mass or its other properties, like 
annihilation cross section or decay branching ratios of the annihilation products. As a bonus, it may be possible to 
infer the actual shape of the DM density profile. 

A related important test of the DM origin would be a measurement of the total 7-ray flux from the GC {^p — 0) 
as a function of the total solid angle Ail. Again, in order to remove the dependence on the WIMP properties, we 
introduce the ratio of the total fluxes from the GC as a function of Afi, 

$^(Ar!- 10-5sr)' 

where the choice of Afi = 10~^sr reflects Fermi LAT's angular resolution but is otherwise fairly arbitrary. 

The quantity is shown in Fig. [31 for our three representative halo models. (The curve for the Klypin, et al., 
model exhibits a non-smooth behavior because of numerical integration over tabulated points, as there is no analytic 



^ We thank G. Johannesson for his comments on this point. 
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expression available for the profile.) Again, since the 7-ray flux from DM annihilation is proportional to p^, it would 
be hard to mimic the behavior with conventional astrophysical sources. In addition, since the ratio falls off more (less) 
rapidly for a more (less) cuspy DM profile, by observing the approximate behavior of the ratio of the measured 
total fluxes one could attempt to infer the cuspiness of the DM halo profile in the area of the GC. The quantities 
introduced in Eqs. ([8]) and ([9]) are not independent as are both proportional, at some level, to the square of DM halo 
density. Nevertheless, they are not the same as the latter involves integration over the photon energy and also over 
the solid angle Afl centered on the GC. 




— Einasto 

— NFW 

,g — Klypin et al 
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FIG. 3: The ratio of the total 7-ray flux ^^{A^l) and ^-,{An = 10~=sr) from DM annihilation in the GC (■!/> = 0) versus 
for three representative halo models. 



IV. CONSTRAINTS ON DM PROPERTIES FROM FERMI LAT MID-LATITUDE DATA 

Next we move on to the discussion of the preliminary data on mid-latitude 7-ray fluxes and the ensuing implications 
for DM mass and distribution. Last Spring, the Fermi Collaboration released preliminary data on mid-lattitude 7-ray 
flux from the region of the sky bounded by 10° < |6| < 20° and < Z < 360° Q (hereafter called Afimid-iattitude for 
brevity), in the photon energy range between 100 MeV and 10 GeV. The spectrum revealed two important features. 
Firstly, it did not confirm the spectrum measured over the same area by EGRET, and is instead softer, with an 
integrated intensity Jlat(^^7 > 1 GeV) = 2.35 ± 0.01 x 10"^cm~^ sec~^ sr~\ compared to Jegret(£'t, > 1 GeV) — 
3.16 ± 0.05 X 10~^cm~^ sec~^ sr~^, where only statistical errors are shown and in both cases the contribution from 
point sources has not been removed. 

The second important feature of the Fermi LAT spectrum is that it approximately agrees with the estimate of the 
contributions astrophysical contributions, as computed by GALPROP (compare the right panel of Fig. 1 of Ref. ^) 
which, in considering DM signatures, we will treat as background. 

It is clear that the good agreement between the measured flux and the astrophysical background puts some con- 
straints on the allowed contribution from DM annihilation. This is illustrated in Fig. |4] where we plot the quantity 
E'^d^j/dEj from DM pair annihilation computed using DarkSusy [T^ for a neutralino with mass — 25 GeV (solid), 
50 GeV (dash) and 100 GeV (dot-dash) and the Klypin, et al. (blue) or the Einasto (black) profiles, and averaged 
over Arimid-iattitude- In Order to reproduce av — 3x 10~^^cm'^/s, typical of thermal WIMPs, we chose the underlying 
parameters of the general Minimal Supersymmetric Standard Model (MSSM) as follows: /i — —M2 — 68.4 GeV, 
ruA = 116.1 GeV for = 25 GeV, 11 = -M2 = 113.3 GeV and tjia = 166.2 GeV for = 50 GeV, and 
II = 211.9 GeV, M2 = -205.9 GeV and tua = 144.2 GeV for = 100 GeV, in addition to fixing tan/3 = 10, 
the average squark mass fhq — 2.5 TeV, At = 1.13 TeV and Ai, = —1.92 TeV. For comparison, we show the Fermi 
LAT data. It is clear that, even for rather small and rather cuspy profiles, the DM contribution is at least an 
order of magnitude too small. 
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FIG. 4: The quantity of 7-ray flux from DM annihilation only averaged over 10° < |&| < 20° and < Z < 360° versus Ej 

for — 25 GeV (solid), 50 GeV (dash) and 100 GeV (dot-dash) and the Klypin, et al., model (blue) and the Einasto model 
with best-flt value of a = 0.17 (black). For comparison, the Fermi LAT data points are marked in red, and the associated la 
errors in black. 



On the other hand, for small the DM contribution would be too big for much more cuspy profiles since it scales 
with p^. This allows us to derive an upper bound on such combinations. This is presented in Fig. [5] where on the 
horizontal axis we plot the DM mass and on the vertical axis we plot the upper limit on (J) (mid — lattitudc), which 
is the quantity J{ip) of Eq. ([2]) averaged over 10° < \b\ < 20° and < I < 360°. In deriving the upper limit we first 
establish the amount that DM annihilations can contribute to the Fermi LAT measurement of E'^ averaged over 

the mid-latitude region Ailmid-iattitudc- We do this by reading out from the right panel of Fig. 1 of Ref. [Ej the Fermi 
LAT data and the astrophysical background, and taking their difference. Next, for a given value of the neutralino 
mass we use Eq. ([1]) to compute, up to a normalization constant, (J) and ensure that it does not exceed the 
allowed DM contribution. For definiteness we take the photon energy in the range 0.1 GeV < < 10 GeV, which 
is the same as in the Fermi LAT data. We scan over the MSSM parameters fj,, M2 and tua in looking for the most 
conservative limit, while ensuring that crw ~ 3 x 10~^^cm^/s. The spread in dN^/dE^ is only a few for niy^ > 50 GeV 
but it increases into a few orders of magnitude at lower m^, which explains the rise in the upper limit on (J), shown 
in Fig. O as a solid magenta curve. On the other hand, the DM contribution to the diffuse flux drops down roughly 
as and at large m-^ the upper limit again becomes much weaker. Varying E^^^^ of the Fermi LAT data and 

of the background within their respective error bars moves the upper limit on (J) (mid — lattitude) up and down by 
a factor of a few but it still remains above the values (J) (mid — lattitude) for the Klypin, et al., (blue), the Einasto 
(red), the NFW (black) and the Burkert profile (cyan), which we have plotted for comparison. Clearly, the upper 
limit, at present, puts a rather mild, but non-trivial, constraint on the DM halo density cuspincss. 

It is clear from Fig.[5]that the Fermi LAT data does put a constraint, albeit a rather mild one, on some combinations 
of DM mass and halo profile at mid-latitudes. As more data is accumulated and analyzed from larger areas of the sky, 
especially towards the Galactic Center, and assuming an adequate agreement with astrophysical contributions, the 
upper limit on (J) as a function of the DM particle mass may become stronger. Nevertheless, it is interesting that, 
even with the current data, one can already place some non-trivial constraints on the DM halo profile as a function 
of its mass. Incoming data from the region of the GC is likely to allow one to put a stronger limit on the cuspiness, 
assuming of course no significant deviation from known astrophysical contributions. 
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FIG. 5: An upper limit (solid magenta line) from mid-latitude (10° < |fe| < 20° and < I < 360°) Fermi LAT data on 
(J) (mid — lattitude) versus the neutralino mass taking the photon energy range of 0.1 GeV < E~f < 10 GeV. For compar- 
ison, we show as horizontal lines {J{ip)){inid — lattitude) for the Klypin, et al., model (blue), the Einasto model with best-fit 
value of a = 0.17 (red), the NFW model (black) and the Burkert model (cyan), averaged over the same region of the sky. 



V. SUMMARY 



The hunt for a dark matter signal is in full swing, with several different strategies and experiments reaching promising 
detection sensitivities. The LAT instrument on the Fermi Gamma-Ray Space Telescope has already produced very 
high quality data on diffuse 7-ray emission from intermediate Galactic latitudes and is soon expected to provide 
a very high resolution and precision map of the region of the Galactic Center where the density of dark matter is 
expected to be enhanced. 

In this paper we have argued that, by examining an angular distribution of the upcoming data from the GC, it may 
be possible to see a contribution from DM annihilations for a cuspy enough DM profile, and furthermore to get an 
estimate on the cuspiness, independently of the mass and other properties of the WIMPs constituting the DM. This 
could provide an unambiguous signature of the DM in the GC. 

Next, we showed that the recent Fermi data from mid-latitudes already allows one to put an upper limit on the 
cuspiness of DM halo profile as a function of the WIMP mass assumed to be the lightest neutralino of minimal SUSY. 
We are eagerly awaiting a release of a fuller set of data from the GC and from the whole sky over the photon energies 
extending to much larger values, 100 GeV and beyond, which will hopefully allow to shed more light on the properties 
of the dark matter and its distribution in the Galactic halo. 

While in deriving the upper bound in Fig. [S] we have used the neutralino as a well-motivated case, this choice was 
not essential and similar bounds could likely be derived for other stable WIMP candidates. 
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